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By means of ab-initio calculations we investigate the optical properties of pure a-SiN^ samples, with 
x £ [0.4, 1.8], and samples embedding silicon nanoclusters (NCs) of diameter 0.5<d<1.0 nm. In 
the pure samples the optical absorption gap and the radiative recombination rate vary according to 
the concentration of Si-N bonds. In the presence of NCs the radiative rate of the samples is barely 
affected, indicating that the intense photoluminescence of experimental samples is mostly due to 
the matrix itself rather than to the NCs. Besides, we evidence an important role of Si-N-Si bonds 
at the NC/matrix interface in the observed photoluminescence trend. 



Over the last few years silicon nitride alloys (a-SiN^) 
emerged as possible alternatives to silicon oxides (a-SiO^) 
in silicon-based optoelectronic devices. What makes 
SiN-r appealing is its larger static and optical dielectric 
constant with respect to SiC>2 (8.19 and 4.33 for a-SiaN4, 
4.6 and 2.46 for a-Si02, respectively), which prevents the 
leakage of tunneling currents in nanostructured devices. 
Also very important is its higher stability at large electric 
fields and high temperatures (commercial a-SisNu sam- 
ples melt at about 2175 K compared to about 2000 K 
for a-Si02). Moreover, the reduced band gap of a-SiN^ 
with respect to SiO x (4.7 cV for a-Si3N4 and 8.9 eV for 
a-Si02) requires smaller operating voltages to obtain the 
same electron and hole currents. In analogy to the case 
of a-Si02,[f] the photoluminescence (PL) and electrolu- 
minescence (EL) of a-SiN^ was attributed to the pres- 
ence of Si nanoclusters (NCs). [2, 5-7] However, recent 
experiments [9] indicate, consistently with hypothesis for- 
mulated previously, [8] that a-SiN^ is photoluminescent 
on its own, i.e. even when Si-NCs are not present in the 
sample. Clarifying this point is a prerequisite for further 
developments of this material. In fact, a typical strat- 
egy for tuning the optical properties of materials embed- 
ding NCs is to control the size, shape, density, and other 
structural properties of the nanoparticles. Of course, this 
strategy would be ineffective if the PL of a-SiN^ is not 
due to embedded NCs. Establishing the effect of Si-NCs 
on the optical properties of a-SiN^ samples is also im- 
portant in view of applications of this material in others 
fields, such as light-emitting diodes, [2] lasers, [3] and solar 
cells. [4] 

In this letter we try to shed some light on this point 
by computing the optical gap and the finite temperature 
radiative recombination rate (RR) of pure samples (i.e. 



not containing NCs) at various stoichiometries, and of 
a-SiNz samples containing Si-NCs of diameter 0.5 nm 
< d < 1.0 nm. This will also allow us to estimate the 
effect of the size of the NCs on the optical properties of 
the material. 

Let us start the description of our results from the pure 
a-SiN-r systems. Samples of a-SiN 2 . of stoichiometry 
0.4 < x < 1.8 containing as many as 224 atoms are pre- 
pared following a very careful quenching-from-the-melt 
procedure, requiring long simulation times of about 1 ns. 
The simulations were performed using the classical MD 
package CMPTool [10] and the atomic interactions are 
governed by the Billeter et al. force field. [11] More de- 
tails about the simulation protocol are reported in Ref. 
[12], where it was shown that the adopted procedure pro- 
duces samples that are consistent with the experimental 
ones in terms of concentration of various types of defects 
and bond lengths. Typically, the finite temperature con- 
figurations produced by this approach are highly stressed, 
resulting in an electronic structure that is difficult to ana- 
lyze (e.g. it contains a large number of gap states). This 
is not an artifact of the method as in real conditions 
these configurations are averaged over a suitable statis- 
tical ensemble. However, the calculation of the ensemble 
averaged spectrum of samples of this size is unaffordable. 
We take, therefore, a different approach and compute the 
spectrum over the equilibrium configuration correspond- 
ing to one finite temperature configuration sampled along 
the classical MD simulation. This equilibrium config- 
uration is obtained by optimizing the geometry of the 
systems using the real-space DFT package SIESTA. [13] 
Electronic and optical properties of the relaxed struc- 
tures are then obtained by reciprocal-space DFT calcu- 
lations using the Quantum ESPRESSO package. [14] These 
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latter calculations are performed using norm-conserving 
pscudopotcntials within the local-density approximation 
(LDA). An energy cutoff of 60 Ry on the plane- wave ba- 
sis set is imposed, and a gaussian smearing function of 
0.25 eV is applied on the electronic population in order 
to circumvent convergency issues related to the presence 
of dangling bonds. 

The determination of the gap, not an easy task in disor- 
dered systems, [15-17] is performed using two methods: 
the Tauc method, [18] frequently used with both exper- 
imental and computational data, and the Inverse Par- 
ticipation Ratio (IPR) method. [20] The Tauc approach 
assumes that at the absorption onset holds the relation 
hijjN{hhj) oc (fkt! — Eg) 2 , where N(Hui) is the number 
of interband transitions at a given energy, i.e. the spec- 
tral density, and E g is the optical gap. E g can, there- 
fore, be obtained by the fitting of the above relation. [19] 
The IPR measures the degree of localization of an or- 
bital. IPR = 1 means a maximally dispersed orbital, 
IPR = oo a maximally localized one. Assuming that 
highly localized states (defect states) do not contribute 
to the spectrum, we can estimate the gap by considering 
only those states with an IPR below some predetermined 
threshold. [20] In Fig. la we report the values of the gap, 
calculated using both methods, as a function of the com- 
position of the a-SiNz sample. The gaps reported here 
are in nice agreement with experimental observation on 
hydrogen-free samples, [15, 17] beside a systematic under- 
estimation of the calculated gaps of at most 0.5 eV that 
is most likely due to the well known limitation of DFT 
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FIG. 1: (a) a-SiN x gap vs x curve estimated by using the Tauc 
(squares) and IPR (triangles) methods, (b) Si-Si, Si-N, and 
N-N bond densities, as a function of the stoichiometry of the 
sample. 



in describing excited states. [21] In addition, Both meth- 
ods show that E g (x) increases with x up to the stoichio- 
metric composition and then gently decreases. Following 
Manca,[22] this trend can be explained by noticing that 
at x = 1.33 the density of the more energetic Si-N bonds 
is maximum, while in sub and supra-soichiometric sam- 
ples their density decreases in favor of less energetic Si-Si 
and N-N bonds (see Fig. lb) . The IPR gap shows a sud- 
den increase in the range x — 0.8 — 1.0, which is consistent 
with experimental observations. [16, 17] On the contrary, 
the gap curve obtained by the Tauc method has an un- 
expected convex profile over the entire x range. This is 
because the latter method does not take into account the 
orbitals' localization. Despite the above limitation, the 
Tauc method gives a gap with a correct overall trend, 
consistent with experimental and IPR results. 

Moving to the effect of Si-NCs on the gap, we stud- 
ied samples containing Si-NCs of three diameters: d = 
0.5, 0.8 and 1.0 nm. These samples are prepared by plac- 
ing a spherical Si-NC in a cavity of suitable size created 
in an a-SisN4 sample. The a-SiaN4 sample used here is 
larger than those considered above and contains as many 
as 882 atoms. Such a large sample is necessary in or- 
der to avoid the interaction of the nanoparticle with its 
periodic images and so that the bulk structure of the ma- 
trix is recovered within the boundary of the simulation 
box. The so prepared samples are relaxed by a 50 ps 
long classical MD simulations and then further relaxed 
by ab-initio geometry optimizations. At a variance with 
the pure matrix case, on samples embedding small NCs 
the IPR method cannot be applied. In fact the IPR re- 
moves all the NC confined states that may contribute to 
the optical absorption. At the same time, no low-IPR 
(bulk-like) states are present in the matrix gap for NCs 
of this size. As a result, the E g of systems embedding 
NCs is the same as that of the corresponding pure matri- 
ces. On the contrary, the Tauc method, which is based 
on the fitting of the spectral density over a large energy 
domain, keeps track of the NCs states lying in the va- 
lence and conduction bands. Thus, in the latter case the 
resulting E g (d) curve presents a regular trend with the 
size of the nanoparticle (see Fig. 2). As expected, the E g 
decreases with the diameter of the NC. The modest vari- 
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FIG. 2: Tauc gap as a function of the size of the embedded 
Si-NCs. The gap value of the pure a-SiNi.33 matrix is also 
reported and denoted by d = nm. 
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ation of the gap with the NC size is likely determined by 
the presence of defects, that are known to cause a weak- 
ening of quantum size effects. [23] Also, the convex shape 
of the trend is in contrast with the d~ a trend expected 
for systems governed by quantum confinement (QC).[2] 
By comparing our results with those of Park et al.[2], we 
suggest that the QC-likc behavior is recovered for NCs 
of diameter > 1 nm. 

Let us move to the study of the optical efficiency with the 
stoichiometry of the pure samples and with the size of the 
Si-NCs. The optical efficiency depends on the (thermally 
averaged) RR, that we compute according to the method 
described in Ref. [24]. It is worth stressing that in pres- 
ence of defective samples the RR may critically depends 
on the defect type and concentration, rather than the 
bulk optical properties. As a result, in highly defective 
samples the defect sites may govern the total RR. It is 
therefore important to consider samples obtained follow- 
ing the same procedure in order to have a meaningful 
comparison of the RR. The RR at 300, 600, and 900 K 
for the pure a-SiN^ samples and the a-SiaN4 contain- 
ing Si-NCs are shown in Fig. 3a and 3b, respectively. For 
pure a-SiNz samples, at high T the RR increases with the 
nitrogen content up to £=1.6 and then decreases. This 
is consistent with recent experimental observations. [25] 
The larger recombination times (^20 /zs at least) with re- 
spect to experimental /xs-ns-like PL decay times, [26] are 
mostly referable to the lack of H atoms and of the con- 
sequent hydrogenic passivation in our samples. [23, 27] In 
addition, other approximations used in our calculations, 
such as the LDA, the value of the static dielectric con- 
stant, the lack of phonon interaction and non-radiative 
decay channels, indicate that the calculated RR should 
be considered a lower bound. Also, the amorphous na- 
ture of small NCs is a known quenching factor of their 
optical activity. [26] It is worth stressing that despite the 
systematic underestimation of the computed RR, its de- 
pendence on the matrix stoichiometry and NC size should 
not be affected by the above approximations. 
For samples embedding NCs, the RR do not present a 
clear trend with the size of the nanoparticle. By taking 
the T=600 curve as reference, we observe a RR reduction 
(ei=0.5) or increase ((2=0.8,1.0) of up to a factor 2 with 
respect to the pure matrix. Such variation of the RR 
is probably connected to the emergence of different kind 
and number of defects at the Si/Si3N4 interface, rather 
than to the presence of the NC itself. This explanation 
is also compatible with recent observations that hypoth- 
esize a secondary role of NC in the PL emission spectra 
of samples formed by varying the silicon excess (in our 



case i=1.23 at (2=1.0 nm).[9, 28] 

Finally, the present result supports the Dal Negro et 
al.[26] hypothesis that the PL of a-SiN^ samples is 
due to states associated to Si-N-Si bridge bonds at the 
NC/matrix interface. In fact, the number of such bonds, 
and the associated states, increases with the size of the 
nanoparticle and the PL intensity grows accordingly, at 
least in a limited range of d. 

The small dependency of the PL on the NC size seems to 
indicate that hydrogen-free samples containing NCs are 
less suitable for producing Si-based emitters [9], while the 
possibility of tuning the optical absorption energy might 
be exploited in forthcoming photovoltaic applications. [5] 
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FIG. 3: (a) Thermally averaged RR for pure a-SiN x samples 
as a function of their stoichiometry; (b) percentage variation 
of the RR with respect to the pure a-Si^N^, matrix, as a func- 
tion of the NC diameter d (d = nm denotes the pure matrix). 
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